Probing the classical field approximation - thermodynamics and decaying vortices 
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We review our version of the classical field approximation to the dynamics of a finite temperature Bose gas. 
In the case of a periodic box potential, we investigate the role of the high momentum cut-off, essential in the 
method. In particular, we show that the cut-off going to infinity limit decribes the particle number going to 
infinity with the scattering length going to zero. In this weak interaction limit, the relative population of the 
condensate tends to unity. We also show that the cross-over energy, at which the probability distribution of the 
condensate occupation changes its character, grows with a growing scattering length. In the more physical case 
of the condensate in the harmonic trap we investigate the dissipative dynamics of a vortex. We compare the 
decay time and the velocities of the vortex with the available analytic estimates. 



I. INTRODUCTION 

Experimental realization of Bose-Einstein condensation in 
a dilute gas of alkali atoms ^ stimulated a number of the- 
oretical studies of the quantum degenerate, weakly interact- 
ing many body system. While the lowest energy state of the 
celebrated Gross-Pitaevski equation |Q| describes remarkably 
well the properties of the condensate at zero temperature, the 
dynamics of the Bose gas at non-zero temperatures remains 
a challenge. The most successful description of the tempera- 
ture dependence of the condensate oscillation frequencies and 
their damping rates has been obtained using different versions 
of the two gas models [|^]. While useful, these models assume 
from the very beginning that the system at finite temperatures 
consists of two distinct factions: the condensate and the ther- 
mal cloud. A more fundamental approach would deduce the 
very existence of the condensed part directly from the dynam- 
ics of the many body system. In a series of papers, several 
groups ^, ^, ^, proposed a classical field approximation 
to fulfill this goal. In this approximation, the bosonic field op- 
erator is replaced by a classical field. This dramatic simplifi- 
cation is rooted in quantum electrodynamics, where classical 
electric and magnetic fields arise naturally if the number of 
quanta is large in a given mode. The classical field approx- 
imation is a simple and convenient technique describing the 
condensate in dynamical equilibrium with the thermal cloud, 
at temperatures close to the critical one. We have shown that 
the whole isolated system may be viewed as a single classical 
field undergoing nonlinear dynamics leading to a steady state 
1^, 1^] (see also [||, ^). The condensate is defined as the dom- 
inant term in the spectral decomposition of the time-averaged 
single-particle density matrix. Two cases were discussed: the 
box with periodic boundary conditions [^, |[ 0] and the real- 
istic case of a spherically symmetric harmonic trap ||^]. In our 
procedure it is the observation process and the finite detection 
time that allow for splitting the system into the condensate and 
the thermal cloud. The aim of the present paper is to further 
corroborate the details and the applicability of the classical 
fields approximation. 

In Section ^ we briefly review the approximation based on 



the classical fields. We stress again the need for the time av- 
eraging of the single particle density matrix for unambiguous 
splitting of the system into the condensed and uncondensed 
phases. The classical field approximation is introduced with 
the high momentum cut-off. In Section III we study in some 



detail the role of this cut-off in the model of a cubic box with 
the periodic boundary conditions. We also analyze the depen- 
dence of the cross-over energy on the value of the coupling 
constant. Within our approximation, it is a growing func- 
tion of the scattering length. In Section |v[ for the first time, 
we apply the method to a dynamical dissipative process. We 
phase imprint the vortex on a finite temperature, partially con- 
densed Bose gas and track its decay. We conclude with the 
summary in Section 



II. THE METHOD 

We start with a brief description of our approach (some 
other details can be found also in[^ ^ |^, First, we con- 
sider the simplest case of a system of bosonic atoms of mass 
m interacting via two-body forces. The A^-body Hamiltonian 
written in the second quantization formalism has a form: 



H 



f dV^^'Xr)— 4'(r)-H - f dV4''(r)4''(r)4'(r)4'(r), 
J 2m 2 J 



(1) 

where the first term is the kinetic energy and the second one 
is the energy of two-body interactions. In writing the Hamil- 
tonian (|l]) we made a standard assumption: in the low en- 
ergy limit a two-body interaction potential is approximated 
by a zero-range one. The interaction strength g is deter- 
mined by a single parameter - the s-wave scattering length a.,, 
g - Anh^ajm. 4'(r) is the field operator satisfying equal-time 
bosonic commutation relations: 



[4'(r, f), 'i'Ur', t)] = 6(r - r'). 



(2) 



The field operator can be expanded into any complete set of 
single-particle wave functions. One of these basis sets is of 
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particular importance - the one which corresponds to eigen- 
functions of a single-particle density matrix|g|. This is be- 
cause these functions have direct physical interpretation as 
they define coherent modes of the system. These modes 
are related to, typically performed, single-particle measure- 
ments. Note that a theoretical description of correlated mea- 
surements, which involve a simultaneous detection of a few 
particles, should be based on an appropriate multi-particle 
reduced density matrix. The single-particle mode of macro- 
scopic occupation is a Bose condensate. In the ideal gas case 
eigenfunctions of the single-particle density matrix are obvi- 
ously one particle states of the external potential. The sit- 
uation becomes more complex if particles interact. In gen- 
eral, the eigenmodes are unknown. The only exception is a 
stationary system of particles trapped in a box with periodic 
boundary conditions. Here the symmetry of the problem im- 
poses very strong constrains and eigenfunctions of the single- 
particle density matrix are simply plane waves - even in the 
presence of interactions. The field operator is therefore: 

nr) = y ^e-'^'flk, (3) 

where at is a bosonic operator which annihilates a particle 
of momentum k. A full operator solution of the Heisen- 
berg equation originating from the above Hamiltonian is not 
available and some approximations are necessary. Follow- 
ing Bogoliubov, at zero temperature, when all particles oc- 
cupy a single zero-momentum mode, the corresponding an- 
nihilation operator can be substituted by a c-number ampli- 
tude. Go — > VNaQ. This is possible because the commutator 
[So, Sq] = 1 «; is much smaller than the number of particles 
occupying the zero-momentum state. The remaining terms of 
the expansion (Q) represent quantum corrections which in the 
lowest order of approximation are neglected. Let us note that 
the procedure leads to a substitution of the the full field oper- 
ator by a c-number wave function: 

"Vir) y/N'¥(r). (4) 

The condensate wave function satisfies the Gross-Pitaevskii 
equation which is successfully used for description of low 
temperature behavior of the Bose-Einstein condensate: 

This equation preserves both the energy and the particle 
number, i.e. the normalization of the wave function 1 = 
J d^^r *I'*(r)*I'(r). The Gross-Pitaevskii equation can also de- 
scribe trapped condensates when a term describing the exter- 
nal potential is included. 

Our approach is a simple generalization of the above pro- 
cedure. Let us notice first that at relatively high energies (but 
below the critical one) there exists a number of modes, say M, 
whose occupation significantly exceeds unity, i.e. many dif- 
ferent momentum states, up to a certain k^ax are highly popu- 
lated. Consequently, all corresponding annihilation operators 



in the expansion can be substituted by complex amplitudes 
and the remaining terms may be neglected: 

k„,„ nr 

❖(r) ^ Va^^ J-e-'^ 'ak = VA^*I'(r) (6) 

k 

This approach evidently leads to the same Gross-Pitaevskii 
equation for a wave function describing the "macroscopic 
part" of the whole system, or, more precisely, to its finite lat- 
tice version. This kind of approach is frequently used in quan- 
tum electrodynamics. The only difference from the zero tem- 
perature case is that now this wave function does not corre- 
spond to the minimal energy solution. Instead, it is a relatively 
high energy state. As it has been shown in[|[ |[ ^, regardless 
of the particular choice of an initial state ^^(r, f - 0), differ- 
ent than the Gross-Pitaevskii eigenstate, the system evolves 
towards the same stationary state uniquely determined by the 
energy and the particle number. More precisely, occupations 
of all macroscopically populated eigenmodes - A^l^kP 
fluctuate around their mean values [||. The self consistency 
of the c-number approximation requires that all populations 
have to be large, A^k > 1 ■ To this end, the number of classical 
modes M has to be carefully chosen for each energy and par- 
ticle number[|| ^]. This fact can be checked only a posteriori 
when the stationary populations are obtained. The intimate 
link between the energy, the number of particles and the num- 
ber of modes is the essential ingredient of the whole method. 

There is still one point which requires additional discus- 
sion. In fact, the "quasi-stationary" high energy solution of the 
Gross-Pitaevskii equation represents a pure state of the system 
with all particles in the same single state ^^(r, t). Therefore, it 
seems to be totally unjustified to interpret |flkP as relative oc- 
cupations of different coherent modes of the system. The last 
concept is associated with a mixed state. Detailed inspection 
shows, however, that the relative phases of complex ampli- 
tudes flk vary in time rapidly. On the other hand an observa- 
tion process lasts for a finite period of time, typically of the 
order of hundreds of microseconds. It is sufficiently long for 
a "decoherence" of different eigenmodes: 

1 r'+^'/2 

"k,'^k2 = -7- dra^ (T)flk,(T) ~ |flk,l 5k,,k2- (7) 

Jt-At/2 

In the stationary regime the time averaged quantities do not 
depend on time. On the contrary, a temporal wave function 
of the whole system varies very rapidly. Therefore, what is 
being observed is related to a time averaged single-particle 
density matrix p(ri,r2) = *I'*(ri)*I'(r2). It is the observation 
that reduces a pure state to a mixed one. For particles trapped 
in a periodic box we have: 

pa^=^£e''c--^'iflki'. (8) 

k 

Eq.(^ shows that eigenvalues of the time averaged single- 
particle matrix are normalized occupations of different coher- 
ent modes of the system. We want to stress that the averaging 
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procedure is essential for the correct physical interpretation 
of the high energy solutions of the Gross-Pitaevskii equation. 
Note that instead of averaging over time, spatial avergaing can 
be performed. Their importance can be fully appreciated in a 
more realistic case of harmonically trapped Bose condensates 
when eigenmodes of a single-particle matrix are not known. 
It is just diagonalization of the time averaged density matrix 
which gives simultaneously the eigenmodes and the popula- 
tions. Results of our calculations for a harmonically trapped 
Bose condensate at finite energy are presented in[^. The time 
averaged density has a characteristic bimodal pattern: sharp 
peak of condensed atoms embedded in a broad thermal cloud. 
On the other hand, an instantaneous distribution consists of a 
number of irregular spiky structures [l^]. According to our 
numerical calculations, the average width of a single spike is 
about three times larger than the healing length, the distance 
over which the condensate wave function heals back when 
perturbed locally. The characteristic length scale of the irreg- 
ular structure is slightly larger (but of the same order) than the 
grid spacing necessary for description of the system of a given 
energy and particle number As we have already mentioned, 
the self-consistency of the approximation requires a very pre- 
cise selection of a number of modes M what in the present 
case directly translates into the grid spacing. In our approach 
M becomes an important, physical parameter 



^kmax - ■^I'^kmm 1^ must bc largcr then one. We arbitrarily set 
this value to 5. This way, having a stationary value of a relative 
occupation |flk,„„P, we can determine both and g. It occurs 
that the requirement of self-consistency indicates that in our 
calculations the number of particles grows with M while the 
interaction strength g decreases in order to satisfy the condi- 
tion (||). Only by examining the link between all parameters 
of the classical field approximation can we fully interprete our 
data. 
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III. BOSE-EINSTEIN CONDENSATE IN A PERIODIC BOX 

In this section we explore our approach and present some 
important technical details. We consider here the system 
trapped in a periodic box because in this situation the coher- 
ent eigenmodes are uniquely determined by the symmetry of 
the potential. This greatly simplifies the problem. We nu- 
merically solve the Gross-Pitaevskii equation (|]) on a finite 
grid using the Fast Fourier Transform split-operator method. 
The number of grid points is equal to the number of differ- 
ent momentum states. In order to explore the link between 
the number of macroscopically occupied modes, the number 
of particles and the energy we solve Eq.(^ in a 3D geometry 
varying the number of eigenmodes M only. In the calculations 
we keep the same initial state *P(r, f = 0) which guarantees 
that the energy per particle is constant. Similarly, the value of 
the following product is fixed: 

gN - const (9) 

The results might seem surprising: the larger the number 
of modes M (with other parameters constant) the higher the 
occupation of the condensate mode |aoP — > 1- 

This result requires a more detailed discussion. The total 
energy, the number of modes M and the product of gN are the 
only control parameters in our method. The number of parti- 
cles as well as the coupling strength g do not enter the dy- 
namical equation (||) separately - it depends only on the prod- 
uct ^A^. Nevertheless, the condition of validity of the classical 
field approximation allows to determine the above parame- 
ters quite accurately. An occupation of the highest eigenmode 



FIG. 1: Fraction of particles in the condensate |aol' as a function of 
particle number for gN=const. For comparison we show the curve 
given by the equation |aoP = 1 - 8.2/A'''\ 



In Fig.yJ a relative abundance of the condensate atoms is 
plotted as a function of the number of particles. The conden- 
sate occupation can be quite accurately approximated by the 
following analytic expression: 

l«ol' ^ 1 - (10) 

This behavior is in agreement with a recent finding [[TTl]. Lieb 
and Seiringer studied the properties of a true ground state A^- 
body wavefunction of a system of A^ Bose particles interact- 
ing via two-body forces. They have shown that in the limit 
of A^ — > oo with gN = const the relative population of the 
dominant eigenmode of the corresponding single particle den- 
sity matrix aproaches 1 while the eigenmode tends towards a 
ground state solution of the Gross-Pitaevskii equation. Our 
method shows how this limit is approached. 

In the next series of calculations we estimate the 
interaction-induced shift of the cross-over energy for Bose- 
Einstein condensation. We prefer to use the term "cross-over 
energy" instead of "critical energy" as, strictly speaking, the 
critical energy is defined in the thermodynamic limit only. For 
a finite system (like the one we deal with here) several mea- 
sures of the corresponding "cross-over" energy can be defined 
as there is no sharp phase transition. One of them [|l2| is based 
on the rapid change of the probabiUty distribution of the con- 
densate occupation. 
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We study the occupation of the condensate mode A^o- The 
total energy of the system is changed while and g are kept 
constant. This requires an optimization of the number of clas- 
sical fields M to be taken into account: again, we set the occu- 
pation of the least occupied mode to 5 (A^k„,„ = 5). The above 
constraint does not allow for a continuous change of the en- 
ergy though. When a quasi-stationary state of a given energy 
is reached we trace a temporal population of the condensate 
mode A^o(f) for a long time (about 10^ time steps). This way 
we obtain a probability distribution of the condensate popula- 
tion which is presented in Fig.^ In this figure we show two 
characteristic distributions for different energies of the sys- 
tem. One of them (for the energy per particle E = 43.8^£'i) is 
centered around some non-zero mean value of A^o ~ 5000 (the 
grey histogram), while the other (for the energy E - 49.5fiEi) 
is peaked at A^o = 0. Here £1 = Tp'iljff I {2mL^) is the energy 
of the first excited state. 

0.2 
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FIG. 2: Statistics of the condensate occupation in time. The grey 
histogram presents the case slightly below the cross-over energy (en- 
ergy per particle E = 43.^Tia)), while the transparent histogram corre- 
sponds to the situation above the cross-over energy (energy per par- 
ticle E = 49.5/ici;). Note that the highest bar of the latter histogram 
is off the scale of the figure (its value is 0.87). Here A' = 10^ and 

g = 10-". 



Despite the fact that the values of energy for which both 
histograms are plotted are very similar, the rapid change of 
the character of the statistics evidently signifies the onset of 
Bose-Einstein condensation. The two energies give the lower 
and upper limit on the critical energy for a given interaction 
strength. We cannot determine the critical energy with bet- 
ter accuracy because of the previously discussed link between 
the parameters of our method. In the Fig.^ we show the up- 
per and lower limits for the critical energy as a function of the 
coupling strength g. The critical energy grows with g. The 
same must happen to the critical temperature. We do not esti- 
mate the critical temperature as the method based on fitting the 
thermal distribution to the outer wings of the time averaged 
density profiles introduces a large uncertaintityljsj. There is 
still a controversy in the literature [[l^ about the magnitude 
and the sign of the shift of critical temperature due to interac- 



tions. The accuracy of our calculations, at the present stage, 
does not allow for a precise determination of this quantity. 
Nevertheless, they certainly indicate that the critical energy 
(and temperature) in the case of gas trapped in a box grows 
with interactions. The same conclusion for the classical field 
simulations in a box has been reached recently in . 
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FIG. 3: Estimates of the critical energy of the phase transition based 
on the analysis of pairs of histograms like in Fig.^ Total number of 
particles is A' = 10^. 



IV. VORTEX DYNAMICS AT FINITE TEMPERATURES 

Creation of quantum vortices has provided a long-sought 
clear evidence of supefluidity in Bose-Einstein condensates 
of dilute atomic gases Jl4| , p^ . Following the two pioneering 
experiments, various mechanisms of vortex nucleation have 
been investigated [p^. More recently, experiments have fo- 
cused on peculiar properties of vortex lattices [p^. 

In most experiments vortices have been created in pure, 
essentially zero temperature, condensates. Very little has 
been said about finite temperature properties of these systems. 
Madison et al. [|5]|, operating at temperatures below their de- 
tection limit, have measured finite vortex lifetimes of the order 
of 500-1000 ms. Anderson et al. [ p^ have seen no thermal 
damping on the time scale of 1 s. Finite temperature decay 
of vortex lattices has been also investigated [|9|. Similarly, 
few theoretical works approached the problem. Fedichev and 
Shlyapnikov [ ^o| ] show that it is the interaction of the thermal 
cloud with the vortex that provides a mechanism of energy 
dissipation. Due to that, an off-center vortex spirals out to 
the condensate boundary where it decays through the creation 
of excitations. For typical experimental parameters, they esti- 
mate the vortex liftime to be in the range O.I-IO s. Tempera- 
ture dependence of the critical rotation frequency for the cre- 
ation of vortices has been calculated in [pT||. Bogoliubov as 
well as Hartree-Fock-Bogoliubov-Popov theories have been 
used to study the stability of vortex states at finite temperature 
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pjjj i . The role of filling the vortex core by the thermal cloud 
as a stabilizing factor has been emphasized | 
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FIG. 4: Trajectory of a vortex core. The intensity of lines and dots 
fades away with time. 



We have used the classical field method to study the dy- 
namics of vortices in the presence of a thermal cloud. As an 
example, we consider 10^ ^^Rb atoms in a spherically sym- 
metric trap of frequency w = 27rl00 Hz. First, a suitable phase 
pattern (corresponding to a single quantum of circulation) is 
imprinted on a finite-temperature equilibrium state of an inter- 
acting Bose gas. Experimentally, phase imprinting is realized 
by shining a far off'-resonant light on the condensate through 
a glass plate whose opacity varies continuously with an az- 
imuthal angle. As a result, the condensate atoms acquire a 
phase depending on their location [p3[]. We situate the core of 
the vortex off" the trap center We see that the vortex created in 
this way starts moving around the trap center, slowly spiraling 
out to the border of the condensate. We follow the trajectory 
of the vortex core by time-averaging the gas density over short 
time. A sample trajectory of the vortex is shown in Fig.^. The 
equilibrium state used here for the phase imprinting contained 
about 48% atoms in the condensate (its energy being equal to 
49.46 fuij). The initial displacement of the vortex core from 
the trap center is 2.66d. The Thomas-Fermi radius of a con- 
densate consisting of 48000 ^^Rb atoms in our trap is 5.205t/ 
(d - y/H/niLL) - 1 .Q9fim being an oscillator unit of length) - 
note that the vortex disappears reaching roughly this distance 
from the trap center. 

Fig.|| shows the time dependence of the distance of the vor- 
tex core from the trap center. One can see a slow drift towards 
the condensate boundary superposed on top of small fluctua- 
tions coming from the interaction with the thermal part of the 
system. From Fig.|]one can conclude that the radial velocity 




FIG. 5: Time dependence of the distance of the vortex core from the 
trap center. Time expressed in units of a trap period T. 



of the vortex is roughly constant. Note also the vortex lifetime 
T (here equal to about 37 trap periods). 




FIG. 6: Angular velocity of the vortex line (expressed in units of 
a>) as a function of its distance from the trap center For comparison, 
an estimate of Fetter and Svidzinsky for a pure condensate of 48000 
^'Rb atoms in our trap is shown. 



Fig.|| shows the angular velocity of the vortex core as a 
function of its (varying with time) distance from the trap cen- 
ter One can see that the vortex moves faster and faster, in par- 
ticular in the neighborhood of the condensate boundary. This 
result can be compared with theoretical predictions for pure, 
zero temperature, condensates - however, in this case there is 
no dissipation and off center singly quantized vortices in our 
case would follow closed circular trajectories. For pure con- 
densates the angular velocity of the vortex line diverges at the 
condensate boundary [^^. We also see this kind of behavior 
in our data for a partly condensed Bose gas. To be specific, we 
have compared our data with a variational estimate provided 
by Eq.(87) in p4|] (here f2 = 0, as the trap is nonrotating) tak- 
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ing 48000 ^^Rb atoms. Despite notable differences between 
the two situations (presence or absence of dissipation due to 
the thermal cloud) the agreement is rather good even at the 
quantitative level. 



the misprint in Eq.(12) of - the bracketed factor should 
be raised to the power of (-1/2) and not (1/2) [Eq]). 




FIG. 8: Dependence of vortex lifetime t on its initial displacement 
from the trap center x,„j„. Points come from our simulations while 
the solid line corresponds to Eq.(pl|). 



FIG. 7: 3-dimensional density of the gas cloud with a bent vor- 
tex line. The outer shell is a density isosurface near the condensate 
boundary. A plane with a surface density plot, perpendicular to the 
vortex line and through the trap center, is also shown. 

More features of the vortex structure can be seen in Fig.^ 
In this figure a 3-dimensional density of the whole gas cloud, 
averaged over a short time, is depicted. The outer shell is 
an isosurface of density near the boundary of the condensate. 
The whole vortex line is visible in the upper part of the fig- 
ure: one can see that it is bent at the edge of the condensate 
and it is much thicker there than in the center of the trap. The 
former observation agrees with recent experimental and the- 
oretical results [Q. A plane perpendicular to the vortex line 
and through the trap center is also shown - it contains a sur- 
face density plot with a distinct hole at the vortex line. 

So far we have kept the initial displacement of the vortex 
core fixed and simply studied the features of its trajectory. In 
the following, we look at the dependence of the vortex lifetime 
on its initial position in a trap. Fedichev and Shlyapnikov 
calculate that: 



^/JIT 



T = 



ln(R/x„in) , 



(11) 



where R is the spatial size of the condensate (well approx- 
imated by the Thomas-Fermi radius), x„„„ is the initial dis- 
placement of the vortex line from the trap center, T is temper- 
ature and p„ X Q.lm^^^T^^- jtv' is the mass density of the ther- 
mal cloud (for details of the derivation see [pO[ ; however, note 



In Fig.|| we plot the vortex lifetime as a function of its initial 
position (calculated from our simulations) as well as the val- 
ues given by Eq.(|ll]) (note that temperature of our system, for 
a given condensate fraction, has been estimated with the help 
of the ideal-gas formula with the finite-size corrections taken 
into account). One can see that it does agree quite well with 
the analytical estimate. However, one should note that one of 
the assumptions leading to the derivation of Eq.(|ll]) is that the 
condition T » /i holds. In our case T Ijj. - 2.55 which may be 
one of the reasons of the slight discrepances between our data 
and the analytically predicted approximate values. Note that 
the total energy of the system is still fixed to 49.46 tio). 

We have also investigated the variation of the vortex life- 
time with the total energy (which amounts to its dependence 
on the condensate fraction). We have seen that this relation is 
much weaker than the one studied above. We have also no- 
ticed that even for a fixed initial displacement (which does not 
specify its position uniquely) the lifetime changes from shot 
to shot - in other words, the complex violent behavior of the 
thermal cloud in different regions of the trap makes the vortex 
lifetime a stochastic parameter 



V. CONCLUSIONS 

We have reviewed our version of the classical field approx- 
imation to the dynamics of a finite temperature Bose gas. It 
provides an unambiguous splitting of the system into the con- 
densed and uncondensed parts with the help of the time av- 
eraging of the single particle density matrix. In the case of 
a periodic box potential, we have investigated the role of the 
high momentum cut-ofF, essential in the method. In particular, 
we have shown that the cut-ofF going to infinity limit decribes 
the particle number going to infinity with the scattering length 
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going to zero. In this weak interaction limit, the relative pop- 
ulation of the condensate tends to unity. We have also shown 
that the cross-over energy, at which the probability distribu- 
tion of the condensate occupation changes its character - one 
of several measures of the energies signifying the phase tran- 
sition for the finite systems - grows with a growing scattering 
length. In the more physical case of the condensate in the har- 
monic trap we have investigated the dissipative dynamics of a 
vortex. We have compared the decay time and the velocities 



of the vortex with the available analytic estimates. 

The next step in going beyond the approximations em- 
ployed in our model is the estimation of the influence of quan- 
tum corrections to the classical fields, in particular the study 
of the corresponding time scales. 
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